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In terms of the Painleve-Gullstrand-Lemaitre coordinates 
a rather general scenario for the gravitational collapse of an 
object and the subsequent formation of a horizon is described 
by a manifestly C°°-metric. For a 1+1 dimensional model of 
the collapse the leading contributions to the Bogoliubov coef- 
ficients are calculated explicitely and the Hawking tempera- 
ture is recovered. But depending on the particular dynamics 
of the collapse the final state represents either evaporation 
or anti-evaporation. The generalization of the calculation to 
3+1 dimensions is outlined and possible implications are ad- 
dressed. 

PACS-numbers: 04.70.Dy, 04.70.-s, 04.62. +v. 



I. INTRODUCTION 

One of the great challenges of theoretical physics is 
the quest for an underlying law that unifies quantum 
theory and general relativity. The investigation of quan- 
tum fields in curved space-times is expected to provide 
a chance of achieving some progress towards this aim. 
Fulling's discovery 0] of the non-uniqueness of the parti- 
cle interpretation in curved space-times may be regarded 
as a basis for various fundamental effects, see e.g. |j] |27j] 
and references therein. Perhaps the most prominent ex- 
ample is the Hawking effect M which predicts the evap- 
oration of black holes. There are two alternatives for 
the investigation of this striking effect: Originally Hawk- 
ing calculated the Bogoliubov coefficients via the geo- 
metric optics approximation (backwards ray-tracing) in 
Rcf. jfl. In contrast to this dynamical treatment Un- 
ruh jq] imposed boundary conditions on the state in the 
static regime in order to reproduce the main features of 
the Hawking effect. Since the state defined in this way - 
the Unruh state - is completely stationary, it merely de- 
scribes the late-time part of the radiation. Of course, in 
general there exists some temporary amount of created 
particles that depends on the dynamics of the collapse. 
But according to Ref. (J] the number of these particles 
is finite with the result that they disperse after a finite 
period of time and thus do not affect the (divergent) late- 
time radiation. 

Ergo it appears quite natural to assume that the state 
after the complete gravitational collapse coincides up to a 
finite number of particles with the Unruh state describing 
evaporation - independently of the particular dynamics 
of the collapse. The question of whether this assertion 



is strictly correct will be subject of the present article. 
For that purpose we shall calculate the number of cre- 
ated particles explicitely without employing the geomet- 
ric optics approximation. It will turn out that the above 
statement is justified only for one particular branch of 
the rather general class of dynamics of the collapse. 

This paper is organized as follows: In Section || we 
set up the basic properties of the quantum field un- 
der consideration. A brief introduction into the con- 



cept of Hadamard states is presented in Sec. II A. The 



number of created particles is calculated in Section III 



In Sees. [II A and III B we deduce the eigen modes in 
terms of the Schwarzschild and the Painlcvc-Gullstrand- 
Lemaitre coordinates, respectively. The Bogo liub ov co- 
efficients a rc der ived explicitely in Sees. HI C| and [II D. 
In Section IIIE| the relevant expectation values of the 
energy-momentum tensor are calculated. We shall close 
with a summary, some conclusions, a discussion, and an 
outlook. 

Throughout this article natural units with G = Ti — 
c = fee = 1 will be used. Lowercase Greek indices 
such as /i, v vary from (time) to 3 (space) and describe 
space-time components (Einstein sum convention). Up- 
percase Roman indices /, J denote complete sets of quan- 
tum numbers. 



II. GENERAL FORMALISM 

We consider a minimally coupled, massless and neutral 
(i.e. Hermitian) scalar (spin-zero) quantum field $ prop- 
agating on a globally hyperbolic space-time {M,g^). 
Global hyperbolicity demands strong causality and com- 
pleteness, cf. ^8|. (Without these requirements the time- 
evolution of the quantum system is not well-defined and 
unitary.) In the Heisenberg representation the kinemat- 
ics of the field $ is governed by the Klein-Fock-Gordon 
equation 



□ $ = 



1 



■ du 



-gg^dM =0 



(1) 



Strictly speaking, the quantum field is represented by 
an operator-valued distribution $ and hence the above 
equation has to be understood in this sense: OF = — > 
&[F] = 0. In a globally hyperbolic space-time the wave 
equation (|l|) possesses unique advanced and retarded 
Green functions A^ v (x,x') and A re t(x,x'), respectively. 
Employing these distributions one may accomplish the 
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canonical quantization procedure via imposing the co- 
variant commutation relations 



$>(x),$(x') = A let (x,x') - A adv (x,x') . 



(2) 



The solutions of the equation of motion (|]) obey a sym- 
plectic structure induced by the inner product 



(F\G)=i I F*d^G, 



(3) 



G = Fd^G - Gd^F. 



With the aid of 

Gauss' law one can show that the inner product (j^) 
is independent of the particular Cauchy surface £ for 
any two solutions of the Klein-Fock-Gordon equation 
UF = OG = 0, cf. [§§. It should be mentioned here 
that the measure dT,^ used above already contains vol- 
ume factors like s/—ffs and is normalized according to 
g d^x. 

The canonical commutation relations (0) imply 



d£ M dx^ 



)■( 



= (F\G) 



(4) 



As a result the inner product of the field $ with positive 
(Fj) and negative (Fj) frequency solutions of the Klein- 
Fock-Gordon equation, respectively, with (Fj\Fj) = 
— (Fj\F* T ) = 5(1, J) and = defines creation 

and annihilation operators, respectively. As it is well- 
known, these operators and thus also the associated num- 
ber operators depend on the particular choice of the solu- 
tions -Fj. This ambiguity represents the non- uniqueness 
of the particle interpretation (see e.g. jj]) and may be 
regarded as the basis of the phenomenon of particle cre- 
ation induced by the gravitational field. 

Averaging the operator- valued distributions &(x) with 
n-point test functions B n £ C^(Ai n ) via 

$ n [B n }= Jd 4 Xl --- J d 4 x n $(x 1 )---$(x n ) 

Xfinfe, ■ • ■ ,X n ) (5) 

we acquire well-defined operators $"[B„]. The complete 
set of all these operators (constructed for all test func- 
tions) generates the *-algebra containing all possible ob- 
servables of the quantum system (with the unit element 
1 = 

The states g of the quantum system can be introduced 
as linear g(fJ,X + vY) = fig(X) + vg(Y) and non-negative 
g(Z^ Z) > functionals over the *- algebra with unit norm 
g(l) = 1. All these states g build up a convex set, i.e., 
for any two states g\ and gi also the convex combination 
Q\ = A^?! + (1 — A) 02 with < A < 1 represents an allowed 
state. The extremal points of this convex set correspond 
to the pure states g = \&) Since every convex set is 
the convex hull of its extremal points all (mixed) states 
can be written as a (possibly infinite) linear combination 
of pure states. 



In order to decide whether a state is pure or mixed in 
character one has to consider the complete algebra. Fo- 
cusing on a sub-algebra a pure state may display proper- 
ties that are usually connected with mixed states. This 
observation might be regarded as the basis of the thermo- 
field formalism, see e.g. and ||. 

It might be interesting to illustrate these points 
by some examples: The space-time obeying the 
Schwarzschild geometry possesses a Killing vector medi- 
ating the associated (Schwarzschild) time-evolution. The 
ground state of the quantum field (with respect to that 
Killing vector) in the region outside the horizon is called 
the Boulware j^] state ■ It contains no particles - again 
with respect to the Killing vector measuring the time of 
an outside observer with a large and fixed spatial distance 
to the center of gravity. (A free-falling explorer may well 
detect particles in that state.) Ergo the Boulware state 
is a pure state with respect to the algebra of the exterior 
region g# = I^b) (^bI- The interior domain possesses 



(Again all assertions re- 
re Schwarzschild time.) As 



no ground state at all, cf. 
fer to the Killing field along t 
further interesting states one may introduce the Kubo- 
Martin-Schwinger (KMS, p0|| ) states £>t describing ther- 
mal equilibrium at some given temperature T. Obvi- 
ously these states are mixed in character - at least from 
the exterior point of view. One important example is 
the Israel-Hartle-Hawking ^,|| state gmH which satisfies 
(in the exterior region) the KMS condition correspond- 
ing to the Hawking temperature. It can be shown |l4| ] 
that this state is indeed a pure state with respect to an 
enlarged algebra. The Israel-Hartle-Hawking state gmn 
contains the same number of ingoing and outgoing par- 
ticles (thermal equilibrium). Hence the total energy flux 
vanishes. The phenomenon of the black hole evaporation 
can be described by the Unruh || state g\j. This state is 
defined via two requirements: no ingoing/incoming par- 
ticles/radiation at spatial infinity and thermal outgoing 
radiation near the horizon, see also ]l5| . If one considers 
a gravitational collapse of an object and assumes the ini- 
tial state to be pure in character (e.g. the vacuum) then 
the final state is - of course - also a pure state. (Here 
the notion of a pure state refers in both cases to the com- 
plete algebra.) The question of whether the initial state 
indeed transforms into the Unruh state will be subject of 



Section III 



A. Hadamard states 

In general, the complete convex set is too large and 
contains more states than physical reasonable. One way 
to restrict to physically well-behaving states is to im- 
pose the so-called Hadamard |5l| condition. Hadamard 
states are states for which the symmetric part of the bi- 
distribution 

W [2) {x,x') = g ($(z)$(z')) = Tr j(5$(x)$(x')} , (6) 
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the two-point Wightman [B2| function, obeys the follow- 
ing singularity structure (in a 3+1 dimensional space- 
time) 



l(w^(x,x') + W^ (x',x)) = 
'U(x,x') 



V(x,x!)lns 2 + W(x,x') , (7) 



(27T) 



where symbolizes the principal part. The antisym- 
metric part of W*- 2 - 1 must be consistent with the com- 
mutation relation (||). s denotes the geodesic distance 
ds 2 = g ^dx^ dx v between the space-time points x and 
x' (which is at least in a neighborhood of a regular point 
x unique). The functions U(x,x'), V(x,x') and W(x,x') 
are regular in the coincidence limit x — > x'. Together 
with the normalization U(x, x) — 1 the first two func- 
tions U(x,x') and V(x,x') are uniquely determined by 
the structure of space-time, e.g. V(x,x) = -R^/12 (with 
being the Ricci tensor, see e.g. [pTJ ) . Hence all in- 
formation about the state g enters W(x_, x') only. 

One important advantage of the Hadamard require- 
ment may be illustrated by considering the regularization 
of expectation values of two-field observables, for instance 
the energy-momentum tensor T M „. The Hadamard singu- 
larity structure ensures the validity of the point-splitting 
renormalization technique, cf. |^]. 

It can be shown that for a globally hyperbolic C°° 
(i.e. infinitely differentiable) space-time (M.,guv) the 
Hadamard condition is conserved, i.e., if the two-point 
function has the Hadamard singularity structure in an 
open neighborhood of a Cauchy surface, then it does so 
everywhere [ jjr| . 

If one considers the gravitational collapse of an object 
which can be described by a C°°-metric the above theo- 
rem can be used to deduce the Hadamard condition for 
the final state. (The initial state is assumed to be a reg- 
ular excitation over the ground state and thus satisfies 
the Hadamard requirement. The Minkowski vacuum of 
course meets the Hadamard structure with U = 1, V = 
and W = 0.) On the other hand it can be shown that if 
the state of a field <f> fulfills the Hadamard requirement 
(among further not as strict assumptions, cf. ]l7j]) in the 
whole space-time of the Schwarzschild geometry and es- 
pecially at the horizon then the asymptotic expectation 
values correspond exactly to a thermal radiation with 
the Hawking temperature T = l/(4xR) (see |0, O, 
and [0). Combining the two statements above we are 
able to deduce the Hawking temperature for any collapse 
scenario that can be described by a C°°-metric. 

It might be interesting to discuss the previous consider- 
ations by means of some examples for the Schwarzschild 
geometry. Applying the theorems above to the Boulware 
state, i.e. the ground state, it follows immediately that 
this state cannot satisfy the Hadamard requirement - at 
least at the horizon. Indeed this state is singular at the 
horizon - its (point-splitting) renormalized energy den- 
sity diverges there (*b| Tq |* B ) ren I -oo for r J. R. It 



can be shown that the Boulware state as well as every 
KMS state (with an arbitrary temperature) fulfills the 
Hadamard requirement away from the horizon r > R, see 
p3~f . But only the KMS state corresponding to the Hawk- 
ing temperature, i.e. the Israel-Hartle-Hawking state (af- 
ter a suitable extension) meets the Hadamard structure 
at the horizon, see e.g. ]23| , |l8| ] . However, the initial (ap- 
proximately Minkowski) vacuum cannot transform into 
this state during a gravitational collapse of an object, 
cf. p^| . In contrast to the Unruh state the Israel-Hartle- 
Hawking state represents thermal equilibrium also for 
r | oo and the associated amount of particles and energy 
cannot be produced by a collapse, see also Sees. [II D| and 
VI below. 



III. PARTICLE CREATION 

Within the Heisenberg representation the time- 
evolution of the quantum system is governed by the op- 
erators while the states remain unaffected. Hence the 
investigation of the Hawking effect goes along with the 
question: How many (final) Schwarzschild particles con- 
tains the initial state? In general, this number depends 
on the particular initial state and the initial metric as well 
as the dynamics of the metric during the collapse. Ac- 
cording to the considerations in the previous Section we 
assume a C°°-metric throughout. It can be shown that 
the Hawking effect (i.e. the late-time radiation) is indc- 



pendent of the (regular) initial space-time, see Sec. Ill D 
below. Similarly any finite amount of particles being 
present initially does not alter the assertions concerning 
the H awking effect (see the remarks at the end of Section 
HID below). For that reason we assume the initial state 
to coincide with the (initial) vacuum. In this situation 
the number of final particles can be calculated via the 
Bogoliubov /3-coefficients, see e.g. p5| 



Nj ut = <0 in |^ ut |0 in ) 



1/3 



(8) 



In order to calculate these coefficients we have to de- 
rive the structure of the initial modes Ff 1 after the col- 
lapse and to compare them with the out-solutions Fj ut 
by means of the inner product in Eq. (g) . 



A. Schwarzschild metric 

The particle interpretation in quantum field theory is 
based on the selection of an appropriate time-like Killing 
vector. This choice refers to a certain class of associ- 
ated observers whose time-evolution is generated by the 
Killing field. For the flat space-time example, the Killing 
vector mediating the (Minkowski) time translation sym- 
metry accords to a usual beholder at rest whereas special 
Lorentz boosts represent accelerated (Rindler) explorers. 
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Since, in general, different Killing vectors generate dis- 
tinct particle definitions, the Rindler observer does not 
regard the Minkowski vacuum as empty with respect to 
(Rindler) particles. Instead, he/she experiences a ther- 
mal bath, a phenomenon which is called the Unruh effect 

The time-evolution of a beholder at a large and fixed 
spatial distance to the center of gravity is generated by 
the Killing vector corresponding to the Schwarzschild 
time t. The particles that are measured by such an ob- 
server can be described by positive frequency solutions - 
with respect to that time coordinate - of the Klein- Fock- 
Gordon equation. In contrast the evolution parameters of 
further coordinate representations of the Schwarzschild 
geometry (e.g. the Lemaitre metric) represent different 
explorers (e.g. the free-falling one) in general. 

In terms of the Schwarzschild coordinates t, r, 1?, <p the 
3+1 dimensional metric assumes the well-known form 



,/,-> = ( i - _ W - h - _ 

r J \ r 



-r 2 *? 2 - r 2 sin 2 ■& dip 2 . 



dr 2 



(9) 



where R denotes the Schwarzschild radius and describes 
the position of the horizon. 

Strictly speaking, there exist several definitions of a 
horizon, for example the event, the apparent, and the pu- 
tative horizon, cf . |2^] and |27| . The notion of the event 
horizon refers to the global structure of the space-time 
(asymptotical reachability) whereas the apparent horizon 
can be defined by strictly local considerations (trapped 
surfaces). Together with some additional requirements 
(e.g. asymptotical flatness, cf. J2lj) also the putative hori- 
zon represents a local condition: 'time slows to a stop', 
cf. f27|| . In the space-time of the eternal Schwarzschild 
geometry (see Figs, [j] and |^) all these definitions coin- 
cide, but in a dynamical space-time describing the grav- 
itational collapse of an object (see Figs. |^ and ^ this 
coincidence does not hold in general. Within our inves- 
tigations we always refer to a locally defined horizon - 
such as the apparent horizon. 

As it will become more evident later on, the most in- 
teresting region (with respect to the Hawking effect) is 
the vicinity of the horizon. In order to extract the fea- 
tures that are characteristic for this zone we introduce a 
dimcnsionless variable x v i a 



X 



R 



1 



(10) 



This quantity allows for a Taylor expansion in the vicin- 
ity of the horizon. As another useful tool we define the 
Regge- Wheeler tortoise coordinate 



dr 



1 - R/i 



= R\n X + 0[ X } 



(11) 



For reasons of simplicity we restrict our further consider- 
ations to the (t, r)-sector and drop the angular contribu- 



tions in Eq. (Q). The resulting 1+1 dimensional space- 
time obeys a conformally flat metric when expressed in 
terms of the tortoise coordinate 



ds 2 = 1 



R 



(dt 2 -drt) 



(12) 



As a result the equation of motion (Q) simplifies to 



\dt 2 drl) 



(13) 



in 1+1 dimensions. (In 3+1 dimensions additional term s 
occur and generate slight modifications, see Sec. VII .) 
As a result the positive frequency Schwarzschild 
eigen functions are given by 



-iijjtzniijjr * 



7 i{x)=N° 



x -i&R ii + { x ]) 



(14) 



for r > R and vanish for r < R due to the horizon, cf. [£4j . 
The ingoing and outgoing modes are distinguished by 
£ = ±1. (For t he m odifications occurring in 3+1 di- 
mensions see Sec. VII.) Af out symbolizes a normalization 
factor which may without any loss of generality chosen 
to be independent of £. These eigen functions are rapidly 
oscillating near the horizon which again hints that there 
is the most interesting region: This singular behavior of 
the modes corresponds to the freezing of the kinematics 
of the field (governed by the Klein-Fock-Gordon equa- 
tion) in the vicinity of the (putative) horizon. 



B. Painleve-Gullstrand-Lemaitre metric 

The Schwarzschild metric is quite simple but exhibits 
a coordinate singularity at the horizon and is there- 
fore not C°° there. Hence it is impossible to express 
a manifestly C°° -metric in terms of the Schwarzschild 
coordinates. For this purpose one has to employ other 
coordinate systems. As one possible candidate we 
consider the Painleve-Gullstrand-Lemaitre |?3| coordi- 
nates tpGLj f, c/2. These coordinates emerge from the 
Schwarzschild coordinates ts,r,i!),ip by means of the 
transformation 



dtp GL = dt s + g y R /. T , dr . 

1 — R/r 



(15) 



There exist two branches of these coordinate set distin- 
guished by a = ±1. In the following we shall drop the 
index t — ipQL for convenience. The metric transforms 
into 



ds 2 



1 



R 



dt 2 



■ drdt - dr 



-r 2 d-d 2 



sin 2 $ dip 2 



(16) 
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r=0 



FIG. 1. Penrose diagram of the maximally extended 
Kruskal manifold. Owing to the conformal mapping the null 
lines (light rays) are at ±45°. The horizontal axis of symme- 
try indicates the borderline between future (above) and past 
(below) . Adopting the notation of Hawking and Ellis the fu- 
ture and the past infinity J+ and J— (with t = +00 and 
t — — 00, respectively, as well as r = 00) are denoted by single 
solid lines. Double solid lines symbolize the future and past 
singularity at r = 0. The horizons at r — 2M are represented 
by dotted lines. i° denotes the spatial infinity (with r = 00 
and t finite). £+ and i— symbolize the future and the past, 
respectively, (with t = +00 and t = — 00, respectively, and 
r finite). Representative surfaces of constant PGL time are 
indicated for a = +1. This branch of the coordinates merely 
covers the black hole (future) horizon and singularity - but 
not the white hole (past) horizon or singularity. 




FIG. 2. Penrose diagram of the Kruskal manifold 
with representative surfaces of constant PGL time for the 
branch a = — 1. Obviously this figure can be obtained by 
time-reversing the previous diagram. The cr = — 1 branch of 
the PGL coordinates covers the white hole (past) horizon and 
singularity. 

Although the Painleve-Gullstrand-Lemaitre (PGL) 
metric does not belong to the well-known and frequently 
discussed representations of the Schwarzschild geometry 
(such as the Kruskal, Eddington-Finkelstein, Novikov, or 
Lemaitre coordinates) it possesses several advantages: 

In contrast to the Schwarzschild form the PGL metric 
(and its inverse as well) is C°° except at the singularity 
at r = 0. 



The hyper-surfaces of constant PGL time dt = are 
equivalent to flat Euclidean spaces; the space-time cur- 
vature is encoded in the shift vector (employing the 
Arnowitt-Deser-Misner (ADM) notation). 
The evolution parameter, i.e. the PGL time, corresponds 
to a Killing vector leading to a stationary metric. This 
fact simplifies the particle definition via positive fre- 
quency solutions. 

Asymptotically r f 00 the PGL representation coincides 
with the Minkowski metric (similar to the Schwarzschild 
form). By virtue of Birkhoff's theorem this coincidence 
persists during the dynamical period of the collapse. 
The radial coordinate r directly corresponds to the sur- 
face of the two-sphere {t, r = const} via Airr 2 . 
Last but not least the effective acoustic metric of the 
sonic analogues (37) of the Schwarzschild geometry equals 
- up to a conformal factor - the PGL form, see e.g. [j38| . 

For further discussions of the properties of the PGL 
metric see e.g. Refs. 34-40 and references therein. For 
example Ref. p5j presents a pedagogical presentation of 
the PGL metric as well as its relation to the Eddington- 
Finkelstein coordinates. 

It should be mentioned here that the PGL coordinates 
do not cover the complete fully extended Kruskal man- 
ifold: E.g., depending on the particular branch (i.e. the 
sign of a) the PGL representation covers either the future 
(black hole) event horizon and the future (black hole) 
singularity for a = +1 or the past (white hole) event 
horizon and the past (white hole) singularity for a = — 1 , 
see Figs. [j] and |. However, both branches possess an 
apparent horizon at r = R = 2M. 

Since we shall calculate the inner product in terms 
of the new coordinates we have to transform the 
Schwarzschild eigen functions, i.e. the out-modes. This 
can be done by simply substituting the Schwarzschild 
time via t s = t PGL - aR\nx + 0[x] in Eq. dlj) 



, x i(°-0«tt(l + O[ X }) . (17) 



F? nt {x)=Af° 



(Again we restrict our considerations to 1+1 dimensions.) 
One observes that the modes with £ = a are no longer 
singular (arbitrarily fast oscillating) at the horizon, only 
those with £ = — a still exhibit this property. As it will 
become evident later on, merely the singular modes with 
£ = — a will contribute to the Hawking effect. 

Employing the Painleve-Gullstrand-Lemaitre coordi- 
nates it is possible to write down a manifestly C°°-metric 
modeling a gravitational collapse of an object and the 
subsequent formation of a horizon 



ds 2 = (1 - f 2 {t, r)) dt 2 - 2af(t, r)drdt - dr 2 



(18) 



with / 6 C°°. Initially t } —00 the metric describes 
a regular object with a (relatively) dilute distribution 
of matter and can be approximated (locally) by the 
Minkowski metric f(t J. — 00, r) = f ln (r) <C 1. For rea- 
sons of simplicity we assume the horizon to be formed 
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at t = 0, i.e. f(t > 0,r > R) = / out (r) = y/Rfr. 
(Note, that we did not impose any conditions on the 
structure of / in the interior region, i.e. beyond the hori- 
zon.) Outside the (spherically symmetric) collapsing ob- 
ject the Birkhoff theorem demands a stationary metric 
f{t,r^R) = y/Rfr. _ 

The Jacobi determinant is simply given by s/—g = 1 
and the metric as well as its inverse are smooth g^, g^" S 
C°° . Of course, this assertion holds true only if we omit 
the formation of the singularity at r = 0. But the region 
beyond the horizon is causally separated from the outside 
domain and (as it will turn out later) irrelevant for our 
purposes. 

Considering the sonic analogues of the Schwarzschild 
geometry the function f(t, r) directly corresponds to the 
time-dependent local velocity of the fluid, cf. |37|-|40|]. 

It should be mentioned here that the knowledge of the 
above metric over a finite period of time is not sufficient 
for determining an event horizon - in contrast to the 
eternal (stationary) metric in Eq. ([l6|). The local metric 
above does also not allow for the construction of a Pen- 
rose diagram (see Figs. |l|||) - this requires the extension 
to the complete space-time. Similarly it does not neces- 
sarily contain a space-time singularity. However, one may 
deduce the existence of an apparent horizon at r = R for 
t > 0. 

In a purely 1+1 dimensional consideration the range 
of the coordinate r in Eq. ( [181) might be chosen arbitrar- 
ily. But in order to keep contact to the 3+1 dimensional 
situation it should be specified according to < r < oo. 
The 3+1 dimensional bouncing-off effect (see also Fig. ||) 
at the origin r = can be simulated in 1+1 dimensions 
by an appropriate boundary condition. As already stated 
in Ref . B , Dirichlet or Neumann boundary conditions or 
every linear combination of them are suitable, see also 
the discussions at the end of the next Section. 



i+ 




i— 



FIG. 3. Penrose diagram of the Minkowski space-time. The 
dashed line symbolizes the (regular) origin r = 0. Light rays 
originating from J— bounce off at the origin and propagate to 
J+. The Minkowski time equals the PGL time for f(t, r) = 0. 
Again representative surfaces of constant (PGL) time are in- 
dicated. 



r=0 i+ 




i- 



FIG. 4. Penrose diagram of the collapse to a black hole 
as described by the branch a = +1 of the PGL metric with 
an appropriate function f(t,r). As one can infer from the 
indicated surfaces of constant PGL time, the formation of the 
black hole horizon (dotted line) can be described regularly by 
these coordinates (with a = +1). 




i- 



FIG. 5. Penrose diagram of the collapse to a white hole as 
described by the branch a — —1 of the PGL metric with an 
appropriate function f(t,r). After the formation of the white 
hole horizon (dotted line) no light ray originating from J— can 
reach r = 0. The particular structure of this figure is based on 
the (not necessary) assumption that the singularity at r — 
develops at a finite period of PGL time after the horizon has 
been formed. Again one may infer from the indicated surfaces 
of constant PGL time that the formation of the white hole 
horizon can be described regularly by the branch a = — 1 of 
the PGL metric. In contrast to Figs. 1 and 2 this diagram 
cannot be obtained by time-reversing the previous figure. 

Although the two distinct branches a = ±1 of the sta- 
tionary Painleve-Gullstrand-Lemaitre metric in Eq. ( |l6| ) 
are related to each other via a simple change of the coor- 



G 



dinates or the time-inversion T : t — > —t, the distinction 
between the different collapse dynamics for a = +1 and 
er = — 1, respectively, in Eq. dlq) cannot be removed by 
any transformation. (It is not possible to find a glob- 
ally integrating factor for the differential form.) The two 
branches correspond to two non-equivalent collapse sce- 
narios (see Figs. |] and ||) and - as we shall see later - 
generate completely different final states of the quantum 
field. Nevertheless, in both cases the initial t J. — oo met- 
ric describes a regular object whereas the final t f oo 
metric represents for r > a vacuum solution of Ein- 
stein's equations with a central mass M = R/2. So there 
exists a priori no reason to prefer one of the two branches, 
see also the remarks in Section VI. 



C. Eikonal ansatz 

In order to calculate the Bogoliubov coefficients we 
have to deduce some informations about the in- modes. 
For that reason we adopt the eikonal ansatz and divide 
the field into an amplitude and a phase 

F'^{t,r) = ^=^(t,r)exp{-iwS c (t,r)} 
v ^ 

1" 



x 1 + 



(19) 



Certainly this ansatz will be justified for compact space- 
time domains (which are not too large) with smooth met- 
rics and high (initial) frequencies u> (see also the remarks 
at the end of this Section) . But as it will turn out later, 
this is exactly the limit that is relevant for the Hawk- 
ing effect. Inserting the above expression into the Klein- 
Fock-Gordon equation (|l|) the leading terms in u> govern 
the kinematics of the phase function via 



(d^)g^ (d v St)=0, 



(20) 



This non-linear equation has four separate branches of 
solutions - e.g. for / = one may identify the positive 
and negative frequency solutions on the one hand and 
the ingoing and outgoing components labeled by £ = ±1 
on the other hand 



d t S ( - afdrS^ = £d r S/: , i.e., 
d t St = (af + Z)d r S ( . 



(21) 



However, these four branches will not necessarily be 
separated for arbitrary space-time dependent functions 
f(t, r). E.g., if f(t, r) oscillates with a large elongation, a 
mode which is initially purely ingoing may turn its direc- 
tion into outgoing and so on. Nevertheless, if we assume 
a sufficiently well-behaving dynamics of /, e.g. if it trans- 
forms directly and smoothly from f m to / ou t _ where the 
relevant time-scales are smaller than the length scales (R) 



- the four branches remain separated: In this case the dif- 
ferent branches cannot approach each other close enough 
during the time-evolution. As a limiting case we may 
consider a very rapid change (sudden approximation) of 
the metric f(t,r) » f in (r)Q(-t) + / ou tW©(*)- In this 
situation the final phase function S$ coincides (nearly) 
with its initial form while its time-derivative changes ac- 
cording to Eq. (|2l]) . The sudden approximation does not 
hold in contrast with the high frequency limit since we 
deal with the frequency-independent phase function S^. 

In summary, the above assumption of a rapid col- 
lapse ensures the separation of the four branches, e.g. if 
d r S^ is positive/negative initially then it remains posi- 
tive/negative also after the collapse. The same applies 
to the time-derivative dt St - as long as / < 1, i.e. out- 
side the horizon, see Eq. (pTj). As a consequence the di- 
vision of the modes into ingoing and outgoing - labeled 
by £ = ±1 - can be used throughout. 

Additional complications arise in 3+1 dimensions, but 
the main result - the separation of the four branches - 
persists under appropriate assumptions: 
As demonstrated in Ref. p4j , a regular spherically sym- 
metric 3+1 dimensional space-time without horizon does 
not allow for the definition of ingoing and/or outgoing 
particles. The eigen modes are standing waves, i.e. linear 
combinations of ingoing and outgoing components with 
equal weights. So the bouncing-off effect at r — mixes 
the ingoing and outgoing components during static as 
well as during the dynamical period. In a 1+1 dimen- 
sional consideration this "reflection" may be simulated 
by an effective boundary condition at r = 0, cf. the re- 
marks below Eq. ( |l8| ) . Selecting appropriate coordinates 
the point r = becomes time-dependent. E.g., in terms 
of length and time-scales associated to an outside ob- 
server the center of the collapsing object goes to infinity 
(asymptotically at a null line) owing to the formation 
of the horizon. In terms of these particular coordinates 
the origin r = corresponds to an accelerated mirror. 
Ref. [[10| presents a derivation of the Hawking effect based 
on the moving mirror analogue. 

In contrast, in terms of the Painleve-Gullstrand-Lemaitre 
coordinates the origin r = obeys no time-dependence 
at all. Again we assume the collapse to occur fast 
enough: The metric is presumed to remain stationary un- 
til t = —R (t he beginning of the collapse) and according 
to Sec. [II B the (apparent) horizon at r = R is formed 
at t = 0. Since the Schwarzschild eigen functions vanish 
for r < R it is sufficient to consider the region outside 
the horizon to be formed at r > R. Within this limited 
space-time domain {— R < t < 0, r > R} the ingoing and 
outgoing components are indeed effectively independent: 
It takes every ingoing light ray inside this domain at least 
the time duration At = R (under the assumptions made 
for the function f(t, r) above) to propagate to the origin, 
to bounce off (turning its direction into outgoing) , and to 
reach the radius r = R again. An analogue assumption 
was already imposed in Ref. ||. In such a scenario the 
information about a possible "reflection" at r = cannot 
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influence the relevant region. In summary we arrive at 
the conclusion that (under the assumptions made) the 
four branches in Eq. ( pl| ) are indeed effectively indepen- 
dent also in 3+1 dimensions (see also the discussion at 
the end of Sec. 0). 

It should be mentioned here that the above ansatz 
is not equivalent to the quasi-classical Jeffreys- Wentzel- 
Kramers-Brillouin (JWKB) approximation (expansion 
into powers of ft) - in spite of some similarities. It does 
also not coincide with the geometric optics approxima- 
tion (backwards ray tracing) which was used in Ref. @. 
Instead the eikonal ansatz is based on a consequent ex- 
pansion into inverse powers of the initial frequency lo. 



D. Bogoliubov coefficients 

Now we are in the position to calculate the Bogoliubov 
coefficients and thereby the number of created particles 
explicitely. Unfortunately, it seems to be impossible to 
find a general solution for these overlap coefficients. Nev- 
ertheless, with an expansion into powers of the relative 
distance to the horizon \ and the inverse initial frequency 
1 jio it is possible to extract the leading contribution - the 
Hawking effect. (As it will turn out later, the sub-leading 
parts merely generate finite contributions and thus do 
not affect the late-time radiation.) Per definition the 
Hawking radiation is exactly that part of the radiation 
which persists at arbitrarily late times (if we neglect the 
back-reaction). Hence the number of created particles ac- 
counting for the Hawking effect has to diverge. Any finite 
amount of particles would disperse after a finite period of 
time and cannot generate late-time radiation. (This is a 
consequence of the spectral properties of the wave equa- 
tion. It possesses a purely continuous spectrum and thus 
does only allow for scattering states but no bound states, 
see e.g. Q and |26).) As demonstrated in Ref. Q, the 
divergent number of particles is necessary for the thermal 
behavior in an infinite volume. In order to isolate the di- 
vergent part of the number of created particles we have 
to consider the Bogoliubov /3-coefficients (see e.g. E3]) 



Ft 



3 M F°j ut . 



(22) 



Since the Painlcvc-Gullstrand-Lemaitre coordinates are 
completely regular the measure dT,^ does not contain 
any singularities. As we have observed in the previous 
Sections, the modes Fj n and Fj ut are bounded. In ad- 
dition, the Birkhoff theorem implies that the modes at 
very large spatial distances to the collapsing object are 
not affected by the collapse. Consequently this region 
does not contribute to the /3-coefficients and generates a 
5(co — a/)-term for the a, see also 0]. In summary we 
arrive at the conclusion that all (single) Bogoliubov (3- 
coefficients are finite. As a result the divergence of the 
number of created particles Nj ut must be traced back 



to the summation/integration over the initial quantum 
numbers I — (£, lo) in Eq. (||) 



NT = 



f.\Pu\ 



There are two possibilities for a singularity, the IR- 
and the UV-divergence of the integration over the ini- 
tial frequencies lo. In the limit of small frequencies lo 
the modes become space- and time-independent and ap- 
proach a constant - unaffected by the Klein-Fock-Gordon 
equation. (Here we regard the IR-singular normalization 
as factorized out.) Ergo in the limiting case lo J. 
the in- and out-modes coincide and thus possess a vanish- 
ing overlap with all other modes corresponding to finite 
frequencies. As a consequence the w-integration of the 
(absolute values squared of the) Bogoliubov coefficients 
is IR-save. 

In summary the infinite amount of particles has to 
be caused by the UV-divergence of the integration over 
the initial frequencies in consistency with Ref. [Q. (The 
Hawking effect is dominated by large (initial) frequencies 
only if one considers a fundamental quantum field theory 
without any kind of dispersion. Introducing a cut-off, 
see e.g. Q], as an effective description of some underly- 
ing theory the calculations are different.) 

Recalling the structure of the initial eigen functions in 
Eq. ( |T9| ) we arrive at the conclusion that only singulari- 
ties of the out-modes may induce a UV-divergence. The 
convolution of regular expressions with the for w | oo 
arbitrarily fast oscillating in-modes yields results of or- 
der 1/lo. Ergo the subsequent w-integration would be 
UV-save. Indeed, the out-modes are not regular at the 
horizon - the region that is naturally relevant for the 
Hawking effect. Thus it is adequate to consider the vicin- 
ity of the horizon and the high (initial) frequency limit 
in order to extract the Hawking effect. As it will become 
more evident later, exactly the leading contributions in 
X and 1/lo are sufficient for the derivation of the thermal 
radiation. 

If we choose the Cauchy surface according to £ = {0 < 
r < oo,t — 0} the surface element assumes the form 
dYj^ = (dr, 0) and the /3-coefficients transforms into 



dr F) 



(dt 



-Of 0, 



Fout 
J 



(23) 



with the quantum numbers / = (£,w) and J = (£',a/). 
Inserting the result of the previous Section dtS^ — 
afdr-S^ = we arrive at 



Pu — J dr exp {— icoS^} 



iLO^drSt: -iL0'(l + (jf{cr~- £,'}/ X ) i^-^'R 

V 'lolo' 



xAf(l + 0\x\) 1 + 



(24) 



8 



At this stage the correct meaning of the Landau symbols 
0[x] and 0[l/w] should be explained: All terms which 
are of higher order in \ and of the same order in uj as 
well as all terms which are of lower order in u> and of the 
same order in % - in comparison with the leading con- 
tributions in the integrand above - are neglected. Such 
a detailed consideration is especially necessary for the 
quantity exp{— iluS^} which involves terms like (^w) n . 
These contributions are not neglected - in contrast to 
terms like x(x LU ) n - 

Accordingly, (exploiting the dominance of the vicin- 
ity of the horizon) we may Taylor expand the amplitude 
A tu (t = 0,r) = A tu {t = 0,r = R)+O[x\. The zeroth- 
order term can be absorbed into the overall normalization 
factor Af via 



N — > AfA^ u (t = 0,r = R) , 



(25) 



and the higher order terms are omitted. (Here and in 
the following we do not change the symbol Af for the 
normalization factor and use the same letter also for the 
modified pre-factors.) A similar procedure can be per- 
formed with the phase function S%. But owing to the 
pre-factor lo it is necessary to expand it up to first order 

Ss(r) = S e (r = R) + d r S^(r = R)R X + 0[ X 2 } , (26) 

cf. the remarks after Eq. (|24|). Again the zeroth-order 
term S^(t = 0, r = R) may be absorbed by a redefinition 
of Af 



Af -^Afexp{-iwS^(t = 0,r = R)} 



(27) 



Since we have to integrate over the initial frequency 
u> in Eq. (@) in order to obtain the number of cre- 
ated particles, the remaining unknown first-order term 
d r S^{t = 0, r = R) can be eliminated by a re-scaling of 
the initial frequency 



ti = uj(d r S i (t = 0,r = R). 



(28) 



Of course, such a transformation may be accomplished if 
and only if £,d r S^ is positive. But according to the argu- 
ments at the end of the previous Section the sign of d r S^ 
does not change during the collapse - as long as it oc- 
curs fast enough and regularly. In this situation the four 
different branches of Eq. ( pil| ) do not mix and thus the 
sign of d r S^ equals its initial value, i.e. £. Again we may 
consider a very abrupt change of the metric (sudden ap- 
proximation, see the previous Section) as an illustrative 
example, where the final phase function nearly coincides 
with its initial form. For the Minkowski example it is 
simply determined by d r S^ ~ £ and no redefinition is 
necessary at all. For other initial metrics the redefinition 
of the frequency exactly corresponds to the fact that the 
Hawking effect is independent of the initial (regular and 
stationary) space-time. 

The Jacobi factor arising from the change of the u>- 
integral measure in Eq. (^) again modifies the normal- 
ization Af only. This undetermined normalization factor 



will be fixed later by virtue of the completeness relation 
in Eq. (|34| ) below. After an analogous Taylor expansion 
of the function f(t = 0, r) = 1 + 0[\] we find 



Pu= / d\ exp{-i&Rx} 



oj-uj'a[o--£']/x 



x x l{<J ^' WR ^^ + 0[x])[l + 



(29) 



As expected from the previous considerations, the Bo- 
goliubov /3-coefficients contribute only for a — — £' and 
vanish (in leading order) for a = In that case the 
out-modes are not singular (at the horizon) - only for 
a = — £' they display the arbitrarily fast oscillating be- 
havior. Hence - depending on the sign a - either only 
ingoing (for a = — 1 and thus £' — +1) or only outgoing 
(for a = +1 and thus £' = —I) particles are produced (in 
an infinite amount). 

The integral in Eq. (|29| ) involves generalized 
eigen functions which do not belong to the Hilbert space 
L 2 but are distributions, cf. f2~ij] . Hence it cannot be in- 
terpreted as a well-defined Ricmann integral. But - as 
demonstrated in Ref. Jm| - it is possible to approximate 
(locally) the generalized eigen functions by well-defined 
wave-packets. One way to simulate such an approxima- 
tion is to introduce a convergence factor via x 6 exp{— ex} 
with e | 

For a = the above integral can be solved in terms 
of T- functions. After insertion of the convergence factor 
we can make use of the formula B3| 



/ 



dx ( 



= y- z T(z) 



(30) 




which holds for 5R(y) > and 5R(z) > 0, and - remem- 
bering r(z + 1) = zT(z) - we arrive at 

^ T(2Wi?) (#S>R + £) 2 ^^' R 

(31) 

In view of Eq. ( |30| ) the higher order terms in X ~ i-e. x - 
cause increasing arguments z. Ergo these terms result 
in higher orders in 1/y — i.e. l/u> — consistently with 
our approximation and the arguments at the beginning 
of this Section. In order to evaluate the absolute value 
squared of the /3-coefficient we may utilize the identity 





T(z)T(-z) = - 
to obtain the final result 



z sin nz 



\ft. 



Af_ 

Cj exp{47Tcj'i?} — 1 



1 + 



(32) 



(33) 
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This expression confirms the argumentation at the begin- 
ning of this Section. The remaining integration over lo (or 
Co) is indeed UV-divergent. In addition we observe that 
the terms of higher order in 1/lo (and thus %) that we 
have neglected in our calculations are not UV-divergent 
and hence do not contribute to the Hawking effect. This 
observation provides an a posteriori justification of our 
expansion into powers of 1/lo and x an d the neglect of 
the sub-leading contributions. 

The UV-divergence can be interpreted with the aid of 
the well-known completeness relation 



aua* IK - (3u/3} K = 5(J, K) , 



(34) 



where / again symbolizes the initial quantum number. 
This equality reflects the completeness of the initial 
modes. Special care is required concerning the deriva- 
tion of an analogue expression involving the out-modes 
since some of those solutions are restricted to the region 
inside or outside the horizon, respectively, and these re- 
stricted modes are not complete in the full space-time. 
In order to apply this relation we have to deduce the a- 
coefficients as well. For that purpose we define slightly 
modified Bogoliubov coefficients via 



/8(wV) = 



(35) 



and in analogy the a-coefficient. The modified Bo- 
goliubov coefficients can be analytically continued into 
the complex w'-plane where the relations F* ut (uj') = 
Fout(—w') and hence 



a(u>,u>') — (3(u>,— lo') 



(36) 



hold. This enables us to derive the Bogoliubov a- 
coefficient for lar ge i nitial frequencies lo. Substituting 
lo' — > — u) 1 in Eq. (|3l| ) together with the complex conju- 
gation the only difference between \au\ and \Pu\ is the 
sign in front of the term i^coR. Dividing the absolute val- 
ues of the two coefficients all other terms cancel and the 
convergence factor e determines the side of the branch 
cut of the logarithm in the complex plane. Hence we find 
for large frequencies lo 



= e X p{-27rw / fl} \au\ 1 + O 



(37) 



Inserting Eq. ( ]37| ) into the completeness relation ( pi[ ) and 
considering the(singular) coincidence J = K it follows 



Nj = (0ml iV° ut |0 in ) 



I/* 



IJ 



3 <T,-f 



6.(1,1) 



exp{47rw'i?} - 1 
exp{47ra/i?} - 1 



finite 



finite . 



(38) 



According to the results of Ref. H| the UV-divergence 
of the w-integration of the absolute values squared of the 
/3-coefficients in Eq. ( p3| ) exactly corresponds to the sin- 
gular quantity S^(I,I) — 5~(lo,lo) and thus represents 
the near-horizon (r j R, i.e. j — oo) part AfyV of 
the infinite volume divergence AfyV — AfyV + AfyV of 
the continuum normalization. As explained in Ref. pi] , 
the infinitely large amount of particles is necessary for 
(quasi) thermal behavior in an unbounded volume. 

It is also possible to calculate the Bogoliubov coeffi- 
cients for regular modes (wave packets instead of plane 
waves), cf. Q. In this case no divergences occur and all 
quantities are finite. Thus the (late-time) Hawking ef- 
fect cannot be easily distinguished from the the (collapse- 
dependent or initially present) finite amount of particles 
via isolating the divergent part in this situation. 

As mentioned before, an initial state gi n containing a 
finite number of particles does not change the final results 
concerning the Hawking effect. Inserting the Bogoliubov 
transformation the expectation value counting the num- 
ber of Schwarzschild particles equals the Hawking term 
plus additional contributions in this situation 



Qin [N°j 



N 



Hawking 
J 



(a*jjajK + P*jiI3jk) gin [a\Ak^ 



IK 



(39) 



J A" 



For a state gi n that contains a finite number of initial 
particles the above expectation values vanish in the high 
(initial) frequency limit loi^lok T 00 ■ As a result the / 
and K summations/integrations are not UV-divergent. 
Hence the additional contributions are finite and do not 
affect the (divergent) Hawking effect. E.g., if we as- 
sume the collapsing object to be enclosed by a (arbi- 
trarily large but finite) box with Dirichlet boundary con- 
ditions we may describe an initial thermal equilibrium 
state via the canonical ensemble. In view of the previ- 
ous arguments we arrive at the conclusion that any ini- 
tial temperature does also not affect the final (Hawking) 
temperature in this scenario. (Dropping the assumption 
of a finite box enclosing the collapsing object the situa- 
tion becomes more complicated. In this case the number 
of particles being present initially diverges owing to the 
infinite volume. Hence the usual spatial infinity part of 
the infinite volume divergence still corresponds to the ini- 
tial temperature whereas the near-horizon part obeys the 
Hawking temperature, cf. ]24| .) 

With the aid of similar arguments one can show that 
the Hawking effect - i.e. the late-time radiation - is also 
independent of the initial metric (as long as it is regu- 
lar). The number of particles created during the transi- 
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tion from one to another regular metric is finite. These 
particles disperse after some finite period of time and 
do not affect the (divergent) late-time part of the radi- 
ation in accordance with the arguments in the previous 
paragraph. In terms of the Bogoliubov coefficients this 
degree of freedom exactly corresponds to the redefinition 
of the initial frequency uj. (W e did not need to specify 
the initial metric /i n (0 in Sec. [II B.) 



E. Energy-momentum tensor 

In order to support the conclusions of the previous 
sections concerning the evaporation/anti-evaporation we 
calculate the late-time expectation value of the relevant 
component of the energy-momentum tensor. In con- 
trast to the inherently non-covariant particle concept this 
quantity is manifestly covariant. The general relativistic 
energy-momentum tensor as the source term in the 
Einstein equations can be obtained via variation of the 
action A with respect to the metric. For the scalar field 
under consideration we obtain the well-known expression 



1 fllS 



■JL^- = du§d v §- 9 -^d^cP§. (40) 



The covariant divergence of this tensor vanishes 

V M Tj} = -^=d^ TS) - \ T a ? d v g aP = 0. (41) 

V 9 " 

However, in general this equality does not imply any con- 
servation law due to the exchange of energy and momen- 
tum between the scalar and the gravitational field (sec- 
ond term). Nevertheless, if the space-time possesses a 
Killing vector mediating the time-translation symme- 
try (Noether theorem) we may define a conserved energy 
current 



(42) 



In view of the symmetry of the energy-momentum tensor 
_ rpvn and Eq Q t0 g e ther with the property of 
the Killing vector V M C^ + ^»Cn — this current is indeed 
conserved 



V„ - . 



(43) 



Now we may calculate the energy flux S out of (or into) 
the black (white) hole 



S = J dE» JP = J d£ M (0 in | f |0 in ) Q v , 



(44) 



where X denotes the (cylindrical) hyper-surface enclosing 
the black/white hole. In a 3+1 dimensional space-time 
one may determine X via the Killing vectors mediating 
the spherical symmetry By virtue of the Gauss law the 
above quantity is invariant under deformations of this 



hyper-surface S. Hence we may consider a sphere with 
a radius which is much larger than the Schwarzschild 
radius - where the metric coincides asymptotically (r f 
oo ) with the Minkowski form. In this region the energy 
flux simplifies to 



(45) 



The symmetrization {■, ■} is necessary in order to obtain 
a Hermitian observable T^ v . The minus sign arises from 
gw(r T oo) = — 1. For large radial distances (approxi- 
mately Minkowski) the expansion of the field reads 



r) 



joutjraut^ r )+h.c. 



f 



fi out e -»(H£r) + h c 



(46) 



Insertion of the above expansion into the bilinear form in 
Eq. ([l5]) generates a sum over £ and £' as well as an inte- 
gration over d and u>'. The time-integration in Eq. ( |45| ) 
involves terms such as exp{±icjt ± iu't} and thus gen- 
erates 5(u) ± a/)-distributions. In view of the positivity 
of the frequencies only lo = to' contributes. Similarly 
the remaining spatial dependence exp{±iw(£ — £')r} im- 
plies that merely £ = £' yields relevant contributions 
at large distances r f 00. As a result only one (w,£)- 
summation/integration survives and the late-time radia- 
tion is related to the number of particles via 



-W\ 2 p(0in\Kf |0 ln )u;£. 



(47) 



This relation confirms the conclusions of the previous 
sections: The divergence of (0i„| |0i n ) exactly cor- 
responds to the time-integration and the resulting sin- 
gularity of the 6(u) — w')-distribution. The Bogoliubov 
coefficients and thus also (0; n | N"^ |0i n ) contribute (in 
an infinite amount) only for a = — £. Hence the collapse 
to a black hole described by the branch a = +1 of the 
PGL metric generates an outward (£ = —1) flux at late 
times whereas the collapse to a white hole corresponding 
to a = — 1 leads to an inward (£ = +1) flux at late times. 

IV. SUMMARY 

In terms of the Painleve-Gullstrand-Lemaitre coordi- 
nates it is possible to model a gravitational collapse of 
an object and the subsequent formation of an apparent 
horizon by means of a manifestly C°° -metric. This set of 
coordinates possesses two separate branches (labeled by 
a = ±1). Depending on the particular branch (i.e. the 
sign of a) either only ingoing or only outgoing particles 
are created in an infinite amount. This infinite amount 
of particles obeys a thermal spectrum corresponding to 
the Hawking temperature. 
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V. CONCLUSIONS 



The theorems presented in Section [I A imply that dur- 



ing every collapse scenario that can be described by a 
C°°-metric an infinite number of particles with a ther- 
mal spectrum corresponding to the Hawking tempera- 
ture is created. This statement is verified in the present 
article for a rather general ansatz for a C°°-metric in 
Eq. ([lsl). For that purpose it is neither necessary to im- 
pose any conditions on the metric beyond the horizon nor 
to specify the explicit dynamics of f(t, r) during the col- 
lapse - as lon g as it is regular, i.e. C°°, and fast enough, 

cf. Sec. una 



So the Hawking effect is not the result of a space- 
time singularity but a consequence of the formation of 
a horizon (strictly speaking, an apparent horizon). For 
the derivation of the Bogoliubov coefficients no assertions 
about the metric in the interior region f(t > 0,r < R) 
are necessary at all. In addition, only the modes that 
are affected by the horizon (one-way membrane, cf. |27| ) 
contribute to the late-time (Hawking) radiation, i.e. the 
outgoing particles for the black hole horizon and the in- 
going particles for the withe hole horizon, respectively. 

Thus the properties of the produced particles cru- 
cially depend on the branch of the Painleve-Gullstrand- 
Lemaitre metric under consideration. Adopting the 
Schrodinger representation the two distinct branches gen- 
erate completely different final states g a . Only one 
state represents the phenomenon of evaporation while the 
other state corresponds to anti-evaporation. 



VI. DISCUSSION 

Perhaps the most striking outcome of the presented 
calculation is the fact that - depending on the particular 
branch a of the dynamics during the collapse - the final 
state of the quantum field does not necessarily represent 
evaporation but possibly also anti-evaporation. The phe- 
nomenon of anti-evaporation has already been discussed 
in the literature, see e.g. pUjffi ], but in a different con- 
text (Schwarzschild-de Sitter geometries, see also [pi). 
In contrast the calculation in the present article applies 
to asymptotically flat space-times. 

For one branch the final state coincides - up to a finite 
number of particles - with the 1+1 dimensional analogue 
of the Unruh state q\j describing evaporation. The other 
branch generates the (in some sense) opposite final state - 
corresponding to anti-evaporation. In the following con- 
siderations we shall denote this state as the anti-Unruh 
state p a u for convenience. This state g a \j can be obtained 
from the Unruh state g\j by means of the (Schwarzschild) 
time-inversion T 
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if we regard the (Schwarzschild) metric of the space-time 
as fixed. Induced by the time-inversion T all outgo- 



ing particles turn their direction into ingoing and vice 
versa. Since the neutral scalar field is neither affected by 
the charge conjugation C nor by the parity transforma- 
tion P (in contrast to a pseudo-scalar field) and we con- 
sider a spherically symmetric situation, both, the Unruh 
as well as the anti-Unruh state are not CPT invariant: 
CPTpu — Pali + 1 Qu- These considerations are relevant 
for the investigation of unitarity and time-reversibility, 
see e.g. @@. 

Searching for the physical implementations of the main 
result of the present article there are several possible in- 
terpretations: 

From a conservative point of view one might argue 
that the branch causing anti- evaporation is unphysical 
and should be excluded. This assertion might perhaps 
be supported by physically reasonable constraints on the 
energy-momentum tensor, such as the energy conditions. 
The two branches of the C°°-metric in Eq. dlg| ) - after 
the straightforward generalization to 3+1 dimensions - 
can be used to derive the associated Ricci tensor R^ v . 
Owing to the smoothness of the metric the curvature 
tensor always exists and is C°° as well. By virtue of Ein- 
stein's equations the Ricci tensor reveals the correspond- 
ing energy-momentum tensor which could be compared 
with an appropriate model of a collapsing star or used to 
test the energy conditions, for example. It is well-known 
that appropriate energy conditions exclude the existence 
of some pathological space-time scenarios, such as cer- 
tain worm-holes or time-machines, see e.g. pTj . However, 
one should be aware that the energy conditions may well 
be violated if one incorporates the back-reaction of the 
quantum field. Since the Hawking effect is most relevant 
for small objects and almost negligible at astrophysical 
orders of magnitude one would expect that such quantum 
effects have to be taken into account. 

As another possible interpretation of the result of this 
article one may arrive at the conclusion that black holes 
evaporate but white holes anti- evaporate. The particle 
production by white holes has already been discussed 
in Ref. p2| , but within a different context: The space- 
time under consideration in Ref. [fl2f was obtained via 
the time-inversion of a space-time representing the col- 
lapse of an object to a black hole, i.e. an anti-collapse. In 
contrast the space-time investigated in the present arti- 
cle corresponds to the gravitational collapse of an object, 
cf. Figs. ^ and |5|. Furthermore the initial state in Ref. |l2] ] 
is determined via a factorization assumption which might 
be questioned in general and does definitely not apply to 
the scenario of the present article. As a consequence 
the resulting radiation becomes singular at the retarded 
time of the termination of the horizon - a prediction 
which differs drastically from the outcome of the present 
article. Based on similar scenarios in Refs. Q and jj| 
further instabilities and quantum effects connected with 
white holes are discussed - before (and independently 
of) Hawking's discovery. The comparison of black and 
white holes is potentially interesting in view of the fun- 
damental question of time-reversibility (unitarity and the 
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second law of thermodynamics) of quantum gravity, see 
e.g. fj|. 

The predicted anti-evaporation of white holes is cer- 
tainly relevant for the sonic black/withe hole analogues, 
see e.g. [|37|-|40|] . These flow profiles always possess an (ef- 
fective) black hole and a white hole horizon, see e.g. Fig. 1 
in Ref. |39|. If the fluid accelerates in such a way that 
its local velocity exceeds the speed of the sound (black 
hole horizon) it decelerates below the speed of the sound 
somewhere (white hole horizon) as well. Consequently, 
the derivation presented in this article implies that, if 
the perturbations of the flow profile of the liquid obey 
a quantum field theoretical description (with the result- 
ing quantum fluctuations), then the associated (effective) 
vacuum fluctuations are converted into (quasi) particles 
leading to evaporation for the black hole horizon and to 
anti-evaporation for the white hole horizon. 

In order to discuss the third opportunity regarding the 
interpretation of the outcome of the present article one 
may recall the fact that the Hawking effect is dominated 
by arbitrarily large (initial) frequencies. But at energies 
above the Planck scale one expects the breakdown of the 
treatment of quantum fields propagating in given (ex- 
ternally prescribed) space-times. In the Planck regime 
the back-reaction, for example, should become impor- 
tant. Assertions about the metric in this region (e.g. the 
Planck scale vicinity of the forming horizon) are a very 
delicate issue. Hence one is lead to the assumption that 
the outside observer cannot distinguish the two branches. 
As any explorer at a finite spatial distance to the collaps- 
ing object merely experiences the static Schwarzschild 
metric (Birkhoff theorem), the only possible way to ob- 
tain informations about the collapse is provided by the 
quantum radiation itself. If we now assume that the out- 
side beholder cannot resolve the behavior of the metric in 
the Planck scale vicinity of the forming horizon, he/she 
can obtain no information about the particular branch of 
the metric a priori. Without any knowledge about the 
value of a during the collapse the most natural ansatz for 
the state governing the measurements of an outside ob- 
server is given by - remember the convexity of the states 
discussed in Sec. [n] 

QV + £>aU 



Qo 
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Again we adopt the Schrodinger representation. This 
ansatz complies with the superposition principle of quan- 
tum theory - if one assumes that the back-reaction of the 
quantum fields onto the metric yields relevant contribu- 
tions. 

The above introduced state describes some kind of 
quasi-thermal equilibrium - it contains the same (infi- 
nite) number of ingoing and outgoing particles with a 
thermal spectrum corresponding to the Hawking tem- 
perature. Although the state go displays in 1+1 dimen- 
sions a close similarity to the (1+1 dimensional analogue 
of the) Israel-Hartle-Hawking state £>ihh, in 3+1 dimen- 
sions this quasi-thermal equilibrium state go differs dras- 



tically from the Israel-Hartle-Hawking state, which de- 
scribes (at least with respect to the algebra of observ- 
ables outside the horizon) real thermal equilibrium. The 
expectation value of the number of particles in the Israel- 
Hartle-Hawking state £>ihh exhibits the complete infinite 
volume divergence, i.e. the near-horizon part r* j — oo 
as well as the usual spatial infinity r* f oo, cf. pi| . 
In contrast the analogue expectation value in the states 
go, p y, an d p a u contains the near- horizon part only, see 
Sec. Ill D . As a consequence the renormalized expec- 
tation value of the energy density in the states go, gu, 
and gau decreases for large distances r with 1 jr 2 whereas 
the same quantity approaches a constant value (in view 
of the Stefan-Boltzmann law proportional to T 4 ) in the 
Israel-Hartle-Hawking state £>ihh • 

It might be noted here that - in contrast to the Unruh 
as well as the anti- Unruh state - the state go is CPT in- 
variant: CPTgo = go- Therefore the unitarity and time- 
reversibility problem mentioned above in connection with 
the (anti) Unruh state does not necessarily apply to this 
state. 



In Sec. HI D we have observed that only the region near 



the horizon generates contributions that are relevant with 
respect to the Hawking effect. Exactly the leading terms 
in 1/lo and x gi ye r i se to the UV-divergence account- 
ing for the Hawking effect. The notion of the vicinity of 
the horizon as the region that is essential for the Hawk- 
ing effect may be illustrated via the following gedanken 
experiment: Let us imagine a very thin shell of mat- 
ter with slowly decreasing radius. As long as the radius 
of the shell is larger than the associated Schwarzschild 
radius the number of created and radiated particles re- 
mains finite as a consequence of the regularity of the met- 
ric and the associated eigen modes. If the shell were to 
stop shrinking before it reached its Schwarzschild radius, 
no Hawking effect would be observed. Accordingly, the 
creation of particles accounting for the Hawking effect 
occurs exactly in the space-time region of the formation 
of the horizon. 



In order to support the argumentation in Sec. Ill C 



concerning the independence and separation of the differ- 
ent branches (e.g. cor resp onding to ingoing and outgoing 
components) in Eq. (21) we may consider a conceptual 
clear scenario - where the effective boundary condition 
at r = does not contribute at all - described in the 
following gedanken experiment: At first we suppose a 
small amount of highly charged matter to collapse at the 
center of gravity forming a tiny extreme Reissner black 
hole. The surface gravity of such a black hole vanishes 
with the result that there is no Hawking radiation (at 
this stage). After the formation of the small black hole 
the point r = is hidden by the corresponding horizon. 
Consequently, there is no "reflection" at the origin r = 
in this case. (It is possible to define ingoing and out- 
going particles separately, cf. pi].) If we now suppose 
the matter (enclosing the tiny black hole) to collapse the 
origin cannot generate a mixing of the different branches 
(e.g. ingoing and outgoing). 
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VII. OUTLOOK 

For the calculation of the Bogoliubov coefficients we 
have restricted our consideration to the 1+1 dimensional 
space-time of the (t, r)-sector. Even though the main re- 
sult of this article should persist also in 3+1 dimensions 
(including the angular terms), some additional compli- 
cations arise Q. The Klein- Fock-Gordon equation (jl|) 
assumes for the 3+1 dimensional Schwarzschild metric a 
slightly modified form. After separating the angular de- 
pendence by spherical harmonics the centrifugal barrier 
and curvature scattering effects can be incorporated into 
an effective potential V e s, see e.g. |l7j 

(^-^+Veff(r.,i))&, m = 0. (50) 

V e ff is strictly positive and approaches zero for | 
+oo and for r* J, -oo with 0[l/r 2 } = Q[l/rl] and 
0[x] — 0[exp{r*/R}], respectively. Unfortunately, in 
3+1 dimensions no closed expression (in terms of well- 
known functions) for the eigen modes is available. The 
asymptotic behavior can be derived easily. For r* j. 
— oo the positive frequency solutions again behave as 
exp{— icot ± iwr*} or linear combinations of them. These 
waves are purely ingoing or outgoing, respectively, for 
r* i — oo. But every mode which is purely outgoing 
near the horizon contains for r* | +oo ingoing compo- 
nents as well owing to the scattering at the effective po- 
tential V c ff (inducing transmission and reflection coeffi- 
cients) and vice versa. In order to obtain a complete 
and orthogonal (with respect to the inner product) set 
of positive frequency solutions of the Klein-Fock-Gordon 
equation one has to combine the different opportunities. 
E.g., within the notation of Ref. Q] (see Fig. 1 there) a 
'down' mode is purely ingoing at infinity and (therefore) 
mixed at the horizon. According to the results of this ar- 
ticle one would expect that an infinite amount of particles 
are created by the collapse in this mode - independently 
of the branch of the metric |Q . 

Furthermore the present article considers the most sim- 
ple example of a quantum field theory, i.e. the neutral, 
massless, and minimally coupled scalar field Further 
investigations should be devoted to fields obeying more 
complicated equations of motion. For the spin-zero field 
example one may incorporate potential terms including 
masses M 2 $ 2 or conformal couplings i?^$ 2 /6 and con- 
sider charged (i.e. non-Hermitian) fields. Moreover, it 
would be interesting to extend the examination to fields 
with higher spin, e.g. the electromagnetic field. Never- 
theless, there is no obvious reason why the main conclu- 
sions of this article should not persist. The evaluation of 
the Hawking effect for interacting fields with non-linear 
equations of motion seems to be rather challenging. 

Similarly the space-time under consideration describes 
the most simple example of a black hole. The 
Schwarzschild geometry represents an uncharged and 
non-rotating black hole where the Einstein tensor and 



thereby also the energy-momentum tensor vanish for 
r > 0. The extension of the results presented in this arti- 
cle to more general static (i.e. non-rotating) black-holes - 
e.g. the Reissner solution - seems to be straight-forward, 
see also p4| . In contrast the investigation of rotating 
(i.e. stationary, but not static) black-hole space-times - 
e.g. the Kerr solution - holds more difficulties. 

Apart from the Painleve-Gullstrand-Lemaitre coordi- 
nates there are several other coordinate sets that describe 
the Schwarzschild geometry space-time by a manifestly 
C°°-metric, e.g. the Eddington-Finkelstein coordinates. 
It might be interesting to consider a collapse model in 
terms of these coordinates in analogy to Eq. ( |i~8| ) and to 
compare the results. 

However, one should be aware that all of the previous 
considerations neglect the back-reaction of the quantum 
field onto the metric. So far the quantum field is treated 
as a test field propagating on a given (externally pre- 
scribed) space-time. If one attempts to leave this formal- 
ism several problems arise: The concept of Hadamard 
states as described in Eq. (f7|) is restricted to free fields 
obeying linear equations of motion. The two-point func- 
tion of interacting fields possesses additional singularities 
in general. Consequently - if one regards the treatment 
of quantum fields in classical (general relativistic) space- 
times as a low-energy effective theory of some underlying 
theory - the imposition of the Hadamard condition is 
not obviously justified. Similarly the requirement of a 
smooth C°°-metric may be questioned from this point 
of view. Accordingly, it might be interesting to examine 
the consequences of collapse dynamics that are not C°° 
regarding the Hawking effect [|43| . 

An exhaustive clarification of these problems probably 
requires the knowledge of an underlying law that unifies 
quantum field theory and general relativity. 
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